Noise rectification by a superconducting loop with two weak links 
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We consider a superconducting loop with two weak links that encloses a magnetic flux. The weak 
links are unequal and are treated as Josephson junctions with non- sinusoidal phase dependence. 
We devise a model that takes into account the fluctuation of the critical currents, due to the 
fluctuations of the order parameter in the weak links. These fluctuations are important near the 
onset of superconductivity; in this regime they may significantly weaken and eventually disconnect 
the superconducting loop. As a consequence of these fluctuations and of the resistive noise in the 
junctions, the average dc voltage does not vanish. Our model can be easily extended to provide a 
qualitative description of a recent experiment. 



S3 
O 
O 

Of 

3 ; 



-a 
c 
o 

> 
o 

CO 

in 
o 

CO 

O 

-i— > 
c3 



i 

C 

O 

o 



X 



PACS numbers: 05.40.-a, 74.50.+r, 74.40.+k 



I. INTRODUCTION 

This study is a cross of two well established phenom- 
ena: Brownian motors and fluctuations near a phase 
transition. The offspring of this combination will be a 
circuit that can rectify thermal noise and sustain a dc 
voltage. 

Brownian motors are systems in which thermal noise 
leads to the diffusion of some coordinate (typically, the 
position of a particle) , whereas a time-dependent applied 
potential (typically, periodic in some direction of space) 
confines this coordinate and controls how far and when 
it is allowed to diffuse. If the applied potential is asym- 
metric, then the fluctuating coordinate drifts in a pre- 
ferred direction and Brownian motion is rectified. Many 
reviewsiSi&ai^ and experiments^ on Brownian motors are 
available. 

Brownian motors are also called "flashing ratchets", 
alluding to the fact that the confining potential pulsates 
in time. A related phenomenon is a "rocked ratchet" ; in 
this case the confining potential is static, but an addi- 
tional time-dependent applied potential is superposed to 
it. This addition distorts the confining potential in such 
a way that any of its minima can temporarily disappear 
and the confined object is transferred to the following 
minimum. A rocked ratchet may also operate without 
thermal noise. 

Phase transitions are usually characterized by an "or- 
der parameter" that vanishes for one of the phases and 
has a non zero value for the other phase. Near a sec- 
ond order phase transition, the equilibrium size of the 
order parameter is small and can be comparable to that 
of its thermal fluctuations. The system we will consider 
is a superconducting loop near its transition tempera- 
ture. Having a loop automatically provides the period- 
icity which is encountered in Brownian motors and the 
presence of permanent currents in a preferred direction 
might be thought of as a starting point for the required 
asymmetric potential. 

The feature that in my view turns a superconductor 
into an interesting rectifying system is the presence of 



two competing currents that are governed by indepen- 
dent fields: the supercurrent is governed by the "order 
parameter" , whereas the normal current is governed by 
the electromagnetic field. Both undergo thermal fluc- 
tuations. However, the strength of the electromagnetic 
fluctuations increases with temperature, whereas fluctua- 
tions of the order parameter are most important near the 
critical temperature. This difference gives us the freedom 
to assign independent sizes to each of these fluctuations. 
In the system we will study, electromagnetic fluctuations 
will play the role of "Brownian motion" , whereas fluctua- 
tions of the order parameter will cause the time variation 
of the confining potential. 



The simplest superconducting component we can con- 
sider is a Josephson junction. A Josephson junction is ac- 
tually a zero-dimensional superconducting wire: its state 
is completely described by the values of the order param- 
eter and of the electrochemical potential at its extremes. 
Therefore, in order to prevent the physical ideas from be- 
ing obscured by mathematical complexity, we will model 
our loop by a uniform superconducting wire interrupted 
by Josephson junctions. Modeling weak links by Joseph- 
son junctions will enable us to reduce the dynamics of 
the system to standard textbook procedures, 7,8,9 Ideas 
for using Josephson junctions as rectifying ratchets have 
been studied in Rcfs. Rectification of the mo- 

tion of vortices has been considered in Ref. [Tj. 



In the following section we define the system to be 
considered and the rules that govern its evolution. We 
also provide some motivation for the choices involved in 
the model and some heuristics as to why rectification is 
expected to occur. Section IlIII is a report of the numeric 
results. In Sec. IIVI our model will be compared with 
available experiments. Finally, there is a summary of the 
achieved results. 
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FIG. 1: Schematic representation of the superconducting 
loop. The vertical thick lines represent segments which are 
"strongly" superconducting and the horizontal branches are 
weak links, which are modeled by Josephson junctions. These 
junctions are modeled by "pure" junctions in parallel with a 
resistor and a capacitor. 



II. OUR MODEL 

A. Basic Features 

Figure [I] shows schematically the system we will study. 
1± and 2± are points in the loop. Between 1+ and 2+ 
(resp. — ) there is a strong superconducting segment and 
between 1+ and 1— (resp. 2) there is a weak link. A weak 
link is a portion of a superconducting wire where super- 
conductivity is weaker, typically, a constriction. These 
weak links will be modeled by Josephson junctions. 

The gauge- invariant phase difference across weak link i 
is defined by 7, = — ipt- + (27r/$ ) C A • ^ s and the 

electromotive force by V% = — Hi- + (1/c) J.^ (dA/dt) ■ 
ds, where tp is the phase of the order parameter, $0 
the quantum of flux (contrary to the electron charge, 
$0 > 0), A is the vector potential, fj, the electrochem- 
ical potential, c the speed of light, t the time, and the 
integrals are performed along the link. For junction 1 
the equilibrium value of the critical current will be de- 
noted by ico, the resistance by R\, and the capacitance 
by C\ — C . For junction 2 the equilibrium value of the 



critical current will be denoted by aI c o and, in order to 
reduce the number of parameters in the model, the resis- 
tance will be taken as R2 = R\/a and the capacitance 
as C2 = aC. I c o, Ri and $o/2irc — H/2e will be taken 
as units; accordingly, the units of voltage, time, capac- 
ity and inductance are I c qRi, h/2eI c oRi, hj2el^R\ and 
H/2eI c o, respectively. In these units the ac Josephson 
relation is 



dr/i/dt = Vi 



(1) 



Additional equations are obtained by equating the cur- 
rent around the loop, which is related to the magnetic 
flux encircled by the loop, with the current flowing across 
each junction: 

hi + V, /Ri + ddVi/dt + J Ni = - (- 1) 1 (<p x + 72 - 7l )/L . 

(2) 

Iji is the current that flows across the "pure" junction, 
Vi/Ri the current that ideally flows through the resis- 
tor according to Ohm's law, CidVi/dt is the current that 
goes into the capacitor, and Im is due to thermal fluctua- 
tions; ip x is the magnetic flux induced by external sources 
multiplied by 27r/$ and L is the self inductance of the 
loop. 

Jntj is the usual Johnson-Nyquist current. Let us 
choose a period of time At which is large compared to the 
autocorrelation time, but can be treated as infinitesimal 
in macroscopic processes. Averaged over At, I-pii bas the 
form 



Im = Wi/vAt; 7 N2 = rjg 2 y/a/At , 



(3) 



with 77 

constant, 

random 

average 

({ImAtf 



T 

numbers 
and unit 



2k B T/R ly where k B 
is the temperature, 



is Boltzmann's 
and gi.2 are 
with normal distribution, zero 
variance. (This follows since 
= f At dt> f Q At dt"(I m (t')I Nl (t")) a 



2At f™{I m (Q)I m (s))ds = 2k B TAt/Ri.) 



B. Heuristic Considerations and Additional 
Choices 

If we ignore Ini, Eqs. Q and J2J predict an evolution 
for 7 j that is equivalent to viscous motion down a poten- 
tial (<^ K +72-7i) 2 /2L+Ii+l2, where dli/d-ji = Iji. In 
order to rectify Brownian motion, there should be some 
direction in the 71^-space such that this potential is pe- 
riodic and is not symmetric under reflection. 

This symmetry breaking cannot be achieved if Iji 
and Jj2 are both sinusoidal functions of the phase dif- 
ference, as is the case for tunnel junctions. Fortu- 
nately, for superconducting wires the current-phase de- 
pendence becomes sinusoidal only as a limit (iiiSiili In 
general, Iji has its maximum away from 7, = 7r/2. 
The deviation from sinusoidality can be quite large 
and Iji can even be multivalued. Various types of 
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FIG. 2: Current through a Josephson junction as a function 
of the phase difference for the case considered in Ref.HiL Here 
the unit of current is wc^aH^ /<&o, where w is the cross section 
of the wire, £a the coherence length in the strongly super- 
conducting wire and H c the thermodynamic transition field. 
The length of the weak link is 2d^A = 0.002£a and the curves 
are for several values of F, a parameter that defines the ratio 
between the strengths of superconductivity in both materials. 
The curves are marked with the number log 10 (r/d). Only the 
curve for F = 10 _1 d is nearly sinusoidal. In order to make 
the curve visible, the current values for the case 
were multiplied by 400. 



nonsinusoidal current-phase relation were also encoun- 
tered in junctions such as superconductor-ferromagnet- 
superconductor point contacts^ 

We may gain a better understanding of the present 
model by reviewing Ref. [LJ in some detail. This refer- 
ence considers a long and thin wire A with a region of 
it, denoted by B, replaced by a different material. B 
is the "weak link" . Both materials are almost identi- 
cal, the only difference between them being that B has 
a shorter mean free path. The weak link is characterized 
by two parameters. One of them, denoted by T, is the ra- 
tio between Gorkov's universal function of the impurity 
parameter for both materials. (Denoting the coherence 
lengths by £a,Bj T = {£b/£a) 2 ) The second parameter 
is the length of B, which is written as 2d!; a- The ar- 
ticle provides explicit expressions that exactly solve the 
Ginzburg-Landau equations and permit to evaluate the 
current-phase relation for the wire. We are interested in 
temperature close to critical, so that we should consider 
d <C 1; on the other hand, any value < T < 1 is phys- 
ically admissible. The authors conclude that only in the 
regime d 2 <C V <C d the current-phase relation becomes 
sinusoidal. In order to appreciate the deviation from si- 
nusoidality, Fig. [21 shows the current as a function of the 
phase difference for d = 0.001 and several values of T. 

For simplicity, we shall keep only two harmonics and 
take 

hi (7i) = hi ( sin 7i + Pi sin 2 7») i (4) 
where I C i and /3j are constants that characterize the junc- 



tion i. In order to present results that are clearly visible, 
most of our figures use the values /3i = — 02 = 0.7, but let 
me emphasize that the numerical investigation covered 
the range < |/3j| < 1 and the qualitative conclusions of 
the following section remain valid for small values of |/3j|. 

In Ref. [lOj the second harmonic is brought in by con- 
necting two junctions in scries. It should be noted, 
though, that this method involves the implicit assump- 
tion that the same current flows through both "pure" 
junctions in series; this assumption is not acceptable in 
the context of this article, since the currents through both 
resistors fluctuate independently. 

Using Eq. J3J, the "potential" that describes the evo- 
lution of 7, becomes (<p x + 72 — 7i) 2 /2L — 7 c i(cos7i + 
\fi\ cos 271) — 7 c2 (cos 72 + \02 cos 272). A contour graph 
of this potential is shown in Fig. [3] for a given set of pa- 
rameters. The dark areas are minima, and the light areas 
are barriers between them. If the barriers are high, the 
values of 7^ will be confined within a minimum, but, if 
the barriers are lowered, the noise will cause these values 
to diffuse, either to the right or to the left. Since the 
barrier at the left is closer than that at the right, the 
probability of crossing it is larger than that of crossing 
the barrier at the right. Therefore, if the barriers are 
raised again after a suitable time, the values of 7$ may 
be pushed towards the consecutive minimum at the left, 
with a greater probability than that of being pushed to 
the right. In this way we may expect a non zero average 
drift of 71 + 72 ; 71 — 72 remains always close to tp x . 

We require a mechanism to vary the asymmetric poten- 
tial of Fig. |3 in time. This is provided by fluctuations of 
I c i and I C 2- The maximum current through a junction is 
proportional to \ip\ 2 , where ip is the order parameter. The 
evolution of ip can be described by the time-dependent 
Ginzburg-Landau theory with the addition of a Langevin 
"force" The "standard" terms drive ip towards equilib- 
rium, whereas the Langevin term causes stochastic devi- 
ations from it. Since the Langevin term increases with 
the absolute temperature and the "standard" terms in- 
crease with the distance from the critical temperature, 
the influence of fluctuations is largest close to the critical 
temperature. 

Due to the evolution of tp, the evolution of hi should 
also have the form of a random walk, together with re- 
laxation to the equilibrium value with some characteristic 
time Tj. The explicit implementation we took was 



I cl (t + At) 



n 



Ici(t) + (Sor + S iri )VAt 



Atl c0 



Tl 



At 



(5) 

and, similarly, for I C 2 the index 1 is replaced by 2, I c q 
by aho and 5q by ado- Here At is the short period of 
time we have chosen, ro, r\ and r2 are random num- 
bers distributed between —1 and 1 and the <5j's denote 
the strength of the fluctuations. If hi or I C 2 becomes 
negative, it is reset to 0. 60 describes synchronous fluc- 
tuations of both junctions, whereas S1.2 describe indepen- 
dent fluctuations. If the coherence length is larger than 
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FIG. 3: The values of 7, behave as the coordinates of a particle 
in two dimensions (with anisotropic mass) which feels the 
potential in this graph, in addition to an anisotropic viscosity 
and a Langevin force. In this graph the darker areas represent 
lower potentials. This potential is periodic in 71 +72, with 
period 4ir. The parameters used to evaluate this graph are 
L = Id = Ic2 = l,Pi = -/9a = 0.7, tp x = tt/2. 

the perimeter of the loop, as is the case for mesoscopic 
loops near the critical temperature, we may expect that 
So will dominate, and vice versa. 

If the Si s are very small, Eq. |JSJ reduces to exponential 
relaxation of I C i to its equilibrium value, with time con- 
stant Tf, in the opposite extreme, if the relaxation time 
becomes infinite, we are left with pure random walk. 

In order to reduce the number of parameters in the 
model we took (3\ = — (3% = /3, t\ = r% = r and 62 = oeS%. 
We are still left with quite a large number of parameters, 
but we shall see that several parameters have just a minor 
influence on the results; in addition, several parameters 
can be grouped, in the sense that varying any of them 
has qualitatively the same effect. 

The relaxation times and the size of the random walk 
steps are related by the fluctuation-dissipation theorem. 
For example, let us consider the case that So = and 
the part of the energy of junction 1 that depends on 
I c i has the form k(I c i — I c o) 2 ■ The constants k and 
I c o could be obtained, for instance, from the Ginzburg- 
Landau theory; both depend on the fabrication details of 
the junction. In addition, I c q should be proportional to 
the distance to the critical temperature and n should be 
independent of the temperature. Let F be the "Langevin 
force" that acts on I c \. Then, by standard methods, we 
obtain J x ryo (F(0)F(s))ds = k B T j kt\. From here, it can 
be obtained that the average of the square of the change 
of I c i during the period At will be ksTAt/ kt\. It then 



follows that TiSf = 2/cbT/k. However, since k depends 
on the details of the juction, we will regard r and the Si's 
as independent parameters. 

The precise form of the fluctuations of I C i is not im- 
portant. I have changed the distributions of the random 
numbers r, from uniform to Gaussian, to bimodal and 
to asymmetric (retaining the values of the average and 
the variance) and the difference in the results obtained 
from different distributions is not noticeably larger than 
the typical difference obtained for different runs with the 
same distribution. All that matters is how often and for 
how long the potential barriers in Fig. [3] are low com- 
pared to the noise level. The scattering in the results, 
though, seems to depend on the distribution used. 

One might gain some intuition by comparing the or- 
ders of magnitude of the energies involved in the problem. 
Fluctuations in the electromagnetic field and in the or- 
der parameter are both of the order of fegT; the energy 
required to break superconductivity (~ kIqo) should be 
significantly, but not exceedingly, larger than ksT. In 
this way, only when many random steps accumulate to 
diminish the order parameter, superconductivity is bro- 
ken at a junction. For the parameters used here, I C i 
typically vanished in one out of 10 4 steps. 



III. RESULTS 

The equations above determine the evolution of the 
phase differences 71,2 for given sets of model parameters. 
This evolution was followed, using Euler iterations, dur- 
ing 10 9 steps. The average voltages are then obtained 
by dividing the change of 71^ by the elapsed time. In 
order to cancel out possible biases in the random num- 
bers, we evaluated four sets of phase differences, which 
were obtained by reversing the sign of either /ni or 7n2 
in Eq. J3J), and then took the average of the voltages ob- 
tained for the four sets. The values obtained for V\ and 
V2 were nearly the same, and the reported voltage values 
are their averages. 

Figure 01 shows our main result, obtained for a given 
set of parameters. In spite of the scattering in the results, 
it can be safely concluded that the dc voltage does not 
vanish, that it is a function of the flux, that it vanishes for 
integer and half-integer values of $/<i>o, and has maxima 
in between. 

I have studied the dependence of the voltage on each 
of the parameters of the model. This study has been lim- 
ited to variations of a single parameter each time, with 
the others fixed at their values as in Fig. 0] For Sq, t 
and 77, the voltage is insignificant below some thresh- 
old, then rises to a maximum and decreases slowly. This 
is expected, because for low values of these parameters 
the phase differences have no opportunities to jump be- 
tween consecutive minima of the "potential" , whereas for 
high values there are too many opportunities and ther- 
mal equilibrium is reached. Figure shows the voltage 
as a function of <5q - Qualitatively, the same shape is ob- 
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FIG. 4: Voltage as a function of the magnetic flux ip x /n = 
2<I>/<E>o- (This function is odd, and periodic in $ with period 
$o-) The parameters used here were a = 1, C = 0, L = 1, 
(3 = -0.7, 5 = 0.02, (5i = 0, r = 5 x 10 4 , r\ = 0.45. For the 
dependence of the voltage on the parameters, see text. The 
curve is a guide for the eye. 
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FIG. 5: Voltage as a function of the strength of the fluctua- 
tions of the maximal supercurrent of the junctions. The flux 
is fixed at tp x — 0.5n and all the other parameters (except for 
So) are the same as in Fig.0] The curve is a guide for the eye. 



tained for quantities playing equivalent roles in diverse 
models, e.g. in Fig. 2 of Ref. 19. As a function of T], the 
voltage attains half of its maximum value at 77 ~ 0.3 and 
77 ~ 0.8. As a function of r, the voltage rises sharply to 
a maximum near t = 2 x 10 4 and then decreases very 
slowly. The reason is that the typical distances of I c i 
and I C 2 from their equilibrium values are actually not 
proportional to r, but rather to its square root. 

The dependence of the voltage on the other parameters 
is roughly as follows. It has a maximum close to a ~ 1.1 
and decreases to half value for a ~ 1/3 and a ~ 4. It 
is insensitive to C for C < 1 and decreases to half value 
for C ~ 10 3 . It is also insensitive to L for L < 1 and 
decreases to half value for L ~ 12. For < |/3| < 1 the 



voltage is roughly proportional to — (3 + 0.25/3 3 and, for 
<5i < 0.03, to e" 8051 . 

In principle, the results should not depend on At, but 
in practice the steps cannot be too large if the deriva- 
tives of ji and Vi are approximated by ratios of discrete 
differences; At also cannot be too small, since then too 
many steps are required until I c i,2 become small a sta- 
tistically significant number of times. In our calculations 
we took At = 0.01. In order to estimate the accuracy of 
our results, we repeated some calculations for At = 0.1. 
The voltages obtained were lower by about 20%. 

An important difference between our model and other 
flashing ratchets is that here the variations of I C i in time 
are themselves due to spontaneous fluctuations, and it is 
not obvious who is the agent that invests the work re- 
quired for rectification. Some models have been studied 
in which rectification is performed by fluctuationsjiiLSi 
but those are manifestly nonequilibrium fluctuations, 
whereas our heuristic considerations suggest that in the 
present case rectification is not subject to any detailed 
form of In or of the fluctuations of I C i, so that the ques- 
tion remains open. 



IV. APPLICATION TO EXISTENT 
EXPERIMENT 

Recent measurements on mesoscopic rings composed 
by segments of unequal widths, near the critical 
temperature^ found a dc voltage with apparently the 
same flux-dependence as in our Fig. 01 Thus far we have 
not explained their results, since the voltages we have 
found are smaller than I c qRi by four orders of magni- 
tude, whereas the experimental results are smaller than 
the corresponding product by only one order of magni- 
tude (provided that we identify Ri with the resistance of 
a segment in its normal state and I C i with the maximum 
supercurrent that can flow in it). 

Nevertheless, bearing in mind the differences between 
the experiment and our model, the very existence of a 
non zero flux-dependent dc voltage in our case is already 
remarkable. The experiment differs from our model in 
two essential respects. The first difference is that their 
rings contain only two segments, without the strongly 
superconducting branches (see Fig. Conceivably, by 
assuming that 2jj depends only on jj and by assuming 
that fluctuations affect the size of Iji only (with the shape 
of the functions remaining invariant), essential features 
of the ring behavior are lost. 

The second difference between the experiment and our 
model is that the experiment used an external ac cur- 
rent I x (t) to destroy superconductivity. This feature is 
readily incorporated into our model. Assuming that both 
segments give equal contributions to the self inductance 
of the ring, all we have to do is add I x (t)/2 at the right 
hand side of Eq. (J2J. By adding this term, our system 
becomes a rocked ratchet ^2*2^ 

Figure shows the dc voltages obtained for our model 
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FIG. 6: Average voltage as a function of the amplitude of 
the external ac current. The upper curve was evaluated for a 
period 2-k/uj — 5 X 10 4 , but we also searched for periods in the 
range between 50 and 5 x 10 6 , and no appreciable difference 
was found; the lower curve is for 2n/ui = 10. The other 
parameters used here were a — 2, C = 0, L — 0.03, /3 = 0.3, 
S = 0.02, 5-l=0, t = 800, r\ = 0.45 and tp x = tt/2. 

with an external current of the form I x (t) — I x o sin(wi). 
Some of the parameters used to evaluate these results 
(a, C, L) were estimated from the available experimental 
data; the remaining parameters are arbitrary. The results 
bear close resemblance to those obtained for other peri- 
odically rocked ratchets for rocking periods 2ir/ui 
that are larger by two or more orders of magnitude than 
the characteristic time fi/2e/ c oi?i of the circuit, we ob- 
tain the upper curve, independently of u>. However, for 
periods that are comparable to h/2eI c oRi, we find that 
there is a nearly periodic pattern of rocking strengths for 
which the rectifying effect is present. 

In the experiment, H/2eI c oRx ~ lCU n s, whereas the 
shortest period of the applied ac current was 10 _6 s. The 
relevant curve is therefore the upper one. Indeed, the 
experiment found that the frequency of the ac current 
has no influence. The upper curve in Figgis remarkably 
similar to Fig. 6 in Ref. |2J, and this time the orders of 
magnitude coincide (for the voltages, and also for the 
currents). 

The heuristic argument used in Fig. [3] to predict the 
existence of a dc voltage is not applicable now: in the 
presence of a large external current, the sign of the volt- 
age is opposite to what that argument would predict. 

Even in the presence of an external current, the noise 
and the supercurrent fluctuations play a central role. For 
Ixo = 4 and the other parameters as in the upper curve 
of Fig. H3 setting r\ = Sq = leads to a voltage decrease 
by one order of magnitude; for I x q < 2.1, the voltage 
vanishes when we set r\ = S = 0. 

The experimental temperature is related to the param- 
eters in our model through the values of I c q. The typical 
deviation of 7 C j from its equilibrium value is of the order 
of SiT 1 / 2 ~ {ksT / n) 1 / 2 . A significant response is ob- 
tained when I c q is larger than — but comparable to — this 



typical deviation, i.e., for temperatures slightly below the 
onset of superconductivity, as in this experiment. 

Let us finally discuss two additional experiments in 
which nonuniform superconducting loops lead to a flux 
dependent dc voltage. Long ago de Waele et ali2* mea- 
sured dc voltage on a double point contact SQUID. The 
asymmetry was not in the weak links, as in the present 
article, but rather in the thick parts of the loop: one 
part was a niobium foil and the other a tin needle. No 
controlled driving ac current was supplied to the circuit; 
apparently, the role of the ac current was substituted by 
existing electromagnetic radiation in the lab. When the 
circuit was appropriately shielded, the dc voltage was no 
longer measurable. As in the cases considered here, the 
highest dc voltage was found when the temperature was 
slightly lower than the critical temperature of the weaker 
superconducting material (Sn). The highest dc voltage 
was comparable to the product of the maximal current 
through a point contact times its normal resistance. 

Weiss et alH measured rectification by an YBCO 
SQUID for a wide range of applied frequencies (ac and 
rf). The asymmetry was mainly due to different maxi- 
mal currents of the junctions; this asymmetry leads to a 
rocking ratchet. Their ac results, shown in their Fig. 5(a) 
are remarkably similar to the upper curve in our Fig. 
including orders of magnitude. In later experiments 25 
the same group studied the influence of temperature and 
repeated the experiment for a niobium circuit. However, 
all the temperatures were considerably below the critical 
temperature, so that only variations in our parameter 
r\ were significant, while the fluctuations of I C i had no 
noticeable influence. In addition, they used tunnel junc- 
tions, so that our (3 is also expected to vanish. 

In view of the differences among the reviewed exper- 
imental loops and our "zero-dimensional" model, their 
similar behavior is remarkable and suggests that the fea- 
tures encountered in our results will be present in any 
sort of nonuniform superconducting loop near its transi- 
tion temperature. 



V. SUMMARY 

Superconductors near their transition temperature 
may be regarded as a special class of thermal ratchets, 
in the sense that the fluctuations of its order parameter 
can act as a flashing asymmetric "potential" and thus 
rectify the Johnson noise. We have found a simple model 
that qualitatively reproduces the experimental results of 
Ref. |2lJ For applied frequencies that are comparable to 
the "natural" frequency of the circuit, a complex pattern 
is predicted. In the absence of external current, the volt- 
age decreases by about three orders of magnitude and 
may change sign, but our simulations and our heuristic 
considerations indicate that it does not vanish. 
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